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By Robert L. McKenzie 
Ames Research Center 


SUMMARY 


The chemical kinetics behind normal shock waves in mixtures of carbon 
dioxide and nitrogen are studied. Particular emphasis is placed on a shock 
speed of 8 km/sec in ambient gas compositions containing 10, and 50 percent 
C0 2 by volume. These conditions are typical of those anticipated for entry 
into the Martian atmosphere. Ambient density is varied from 10”6 to 10“ 1 times 
that of the earth sea-level density. The species considered are C0 2 , N 2 , CO, 
NO, CN, 0 2 , C, 0, and N. A system of 17 reactions is initially assumed and 
later reduced to 9 hy removing those found insignificant to the gross thermo- 
chemical behavior. The methods used for estimating exchange reaction rates are 
discussed in detail, and a scheme for correlating experimental activation ener- 
gies is suggested. General features of the nonequilibrium chemistry are pre- 
sented followed by an analysis of the uncertainties due to possible errors in 
estimated reaction rates. Particular consideration is given to the kinetics 
of CN for shock speeds of k to 8 km/sec. The CN concentration is shown to 
rapidly overshoot its equilibrium value for most ambient densities and shock 
velocities anticipated during entry into the Martian atmosphere. 


INTRODUCTION 


Recent studies of high-speed entry into postulated Martian atmospheres 
(e.g., ref. l) have indicated the need for an improved understanding of the 
chemical kinetics in shock heated mixtures of carbon dioxide and nitrogen. 
Numerous investigations have dealt with high temperature nonequilibrium phenom- 
ena in air and its constituents (e.g., refs. 2 and 3)> but relatively few have 
concerned other gases. Howe and Sheaffer (ref. 4) estimated the chemical 
kinetics behind normal shock waves in pure C0 2 , and Bortner (ref. 5) considered 
mixtures of C0 2 and N 2 but made only a cursory estimate of the reaction system. 

This study is intended to provide a more inclusive estimate of the 
expected thermo chemical behavior behind strong shock waves in mixtures of C0 2 
and N 2 based on the eapei-imental reaction-rate data available to date. The 
objective must be limited, however, to establishing only the qualitative 
nature of the post-shock chemical kinetics since current experimental reaction- 
rate data are sparse and many theoretical rate estimates are incorporated. In 
light of the limited experimental data, it is also unjustified to retain any 
consideration of physical phenomena which have less that a major effect on the 
results. Regardless of these unavoidable inaccuracies, the qualitative 



characteristics of these results may still be shown to be -unique. This is 
accomplished by estimating reasonable limits to the errors in the theoretical 
reaction rates and applying them to obtain extreme kinetic profiles. The 
resulting profiles are then compared for all cases. 

The controlling kinetic mechanism of the chemical system is estimated by 
considering all reactions among a prescribed set of species which are thought 
to have possible significance to the gross thermo chemical behavior. Those 
reactions and species which then prove to be of secondary importance are 
eliminated and the small effects of their absence demonstrated. 

The flow-field model chosen is that behind a normal shock wave because of 
its simplicity and direct application to most shock-tube experiments. The 
principal quantities defining the shock-wave conditions are shock speed, 
ambient gas composition, and ambient density. Ambient temperature, whose value 
is incidental to the results, was fixed at 200° K. 1 * A single shock speed of 
8 km/sec was selected for detailed study but lower speeds to k- km/sec are also 
discussed. The speed of 8 km/sec was chosen principally because it represents 
a typical Mars entry condition and also because it serves to define the pre- 
dominant chemical mechanisms for the lower speeds considered. The ambient 
densities considered are 10~ s to 10* -1 times the earth sea-level density. These 
values bracket the high-speed flight regime of aerodynamic interest during any 
planetary entry. The ambient gas compositions selected are 5; 10, and 50 per- 
cent C0 2 (by volume) with the remainder, N 2 . They were chosen in view of the 
three Mars atmosphere models given by reference 6 . 

It is emphasized that the results given by this study are to be considered 
preliminary estimates only. They may be assumed as qualitatively correct, but 
are subject to substantial quantitative adjustment as additional experimental 
results become available . 


SYMBOLS 


A 

D 

E a 

F 

■^t 9 ^r 9 -^v 
% 


frequency factor, defined by equation (8), m 3 /mole-sec 
molecular bond or dissociation energy, J 
endothermic activation energy, J 
partition function per unit volume 

partition functions for translation, rotation, and vibration, 
respectively 

endothermic heat of reaction, J 


1 Ambient temperature affects only the temperature ratio across a strong 

shock wave but has negligible effects on the absolute values of any quantity 

behind the wave front . 
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Planck’s constant, 6 .62 56X10 “ 34 J sec 
equilibrium constant based on molar concentrations 
Boltzmann constant, 1 . 3^05X10 3 j/°K 

reaction rate coefficient in the forward (endothermic) and reverse 
directions , respectively 

total number of significant rotational degrees of freedom 
an arbitrary chemical species acting as a collision partner 
molecular mass, kg/particle 
molecular weight, kg/kg -mole 


Ste/mi, also equal to the ratio: moles of ambient gas mixture per 

mole of local gas mixture 

Avogadro's number, 6. 022X1 0 23 particles/mole 

mole fraction, moles of species j per mole of local mixture 
number of atoms per molecule 

volumetric rate of production of species j by reaction i, 
moles/sec m 3 


species production rate parameter, Q. . 




may 0 \ / P 0 


collision cross section (see eq. (4)) 

total number of vibrational degrees of freedom 

temperature , °K 

ambient temperature , 200° K 

5k 

Ti 

particle time, sec 
laboratory time, sec 
volume 

shock speed, m/sec 

local speed behind the shock wave (shock wave reference frame) 
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V 

y 

y 0 

y 


V 

V s 

distance behind the shock front, m 
reference distance, 1 m 
distance parameter, 


a 

M 


e r 
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M- 


P 

p o 

P 

a 

[A] 


an index representing a vibrational energy mode 
constants in equation (3) 

local concentration of species j, moles of j per mole of ambient 
gas 

characteristic temperature for rotation, °K 

characteristic temperature of the a vibrational mode, °K 

mAm-D 

reduced molecular weight, 

m A -rn B 

mass density, kg/m 3 

reference density, 1.225 kg/m 3 (standard earth sea level) 
p i 

symmetry number in equation (2) 

concentration of species A in moles per unit volume 


Subscripts 

A,B,AB evaluated for the particle A, B, or AB 

e,(e) equilibrium value 

f frozen chemistry value 

i reaction i 

j species j 

(p) peak concentration 

r reactant r 
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o 


reference value 


1 

2 


conditions of the gas 
conditions of the gas 


mixture ahead of the shock wave 
mixture "behind the shock wave 


METHOD OF ANALYSIS 


General Considerations and Assumptions 

This analysis is "based on the following assumptions: A normal shock wave 

propagates at constant velocity into a quiescent mixture of C0 2 and N 2 of -uni- 
form density and temperature. After passage of the shock front the molecular 
internal energy modes of rotation, vibration, and electron excitation reach 
equilibrium with the translational energy prior to the onset of any chemical 
change. The mixture of chemical species mixture just behind the shock, which 
is the same as that ahead, then starts from its so-called "frozen chemistry" 
condition and proceeds through a system of homogenous bimolecular reactions to 
an equilibrium state. All the molecular energy modes follow in equilibrium 
with the local translational temperature. It is well known that some of these 
assumptions, although classical, are not always valid, particularly that of 
locally equilibrated vibrational energies. In view of the qualitative nature 
of this study, however, the addition of any further refinements is unjustified. 

C hemical species .- The complexity of the chemical system to be analyzed 
depends primarily on the number of species considered. The species important 
to the chemical kinetics are assumed in this case to be those significant in 
determining the equilibrium thermodynamic state of the gas. From reference 7; 
these are shown for shock speeds of b to 8 km/sec to be: C0 2 , N 2 , CO, NO, GN, 

0 2 , C, 0, and N. Assuming only these species simplifies this study since, at 
the high temperatures associated with "frozen chemistry" immediately behind the 
shock, the gas tends toward complete ionization. The consequences of this 
assumption should not be serious, however, since at the shock speeds of inter- 
est, the ionization rates are comparable to the chemical rates in an environ- 
ment moving toward negligible ionization (e.g., ref. 8 ). Furthermore, the 
estimates of reference 5 which included ionization reactions and shock speeds 
up to 12 km/sec indicate that even for much higher speeds than 8 km/sec, the 
nonequilibrium concentration of ionized species is not predominant compared to 
the concentration of neutral species. Thus, the consideration of ionized 
species is not expected to significantly affect the conclusions of this study. 
References 9 , 10, and 11 provide the relative abundances of ionized and 
uncharged particles for the frozen and equilibrium thermodynamic states of 
present interest . 

Consistent with the assumption of a locally equilibrated distribution of 
internal energies, specific consideration of the species in internally excited 
states is also omitted. Excited species have been shown to react at higher 
rates than their ground-state counterparts but the effects of such considera- 
tions can be assumed absorbed into the given mean reaction rates which account 
for reactants in all energy states. This is the case when rate coefficients 
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are determined from experiment since most analyses of such, experiments assume 
an equilibrium distribution of energy states except when specific excitation 
reaction rates are given. 

Reactions.- The preliminary set of reactions used for this study is given 
by table I. Each reaction is written with the endothermic direction from left 
to right. The reactions are divided into three general groups, according to 
the manner in which they were studied. The dissociation reactions (No. 1-6) 
are considered fundamental to the reaction system and were included in all 
calculations of the chemical kinetics. In this study they are assumed to 
occur upon collision of the dissociating species with any arbitrary particle, 

M. This assumption is a convenient means of simplifying the reaction system 
although using an arbitrary collision partner rather than a specific one can 
in some cases lead to a factor of 20 error in reaction rate. Because there 
are many other uncertainties which must be incorporated that preclude the pen- 
alties of this one, the benefits of a simplified system make the assumption 
attractive. The air exchange reactions (No. 7-9) are those familiar to most 
nonequilibrium air studies. Their rates have been measured many times and they 
are considered to contribute comparatively little uncertainty to the results 
of this analysis. The remaining reactions (No. 10-17) are those resulting 
from the presence of carbon atoms and molecules containing them. The eight 
carbon reactions given represent less than half the combinations possible from 
the species considered. The others fall into several groups, all of which 
have been omitted. For example, six more reactions may be written which con- 
tain the same reactants as given in table I but result in different products. 

A typical example is shown with reaction 1^, given below with its endothermic 
heat of reaction. 


CO + NO 5: C0 2 + N 102 kj/mole 

The heat of reaction provides an indication of the energy required to activate 
the reaction. These same reactants could also yield 

CO + NO x CW + 0 2 kj/mole 

What appear to be simpler and more likely products in this reaction require a 
much higher energy to occur. The corresponding increase in activation energy 
affects the associated reaction rate in an exponential manner making the reac- 
tion far less effective than its lower energy counterpart. Hence, the higher 
energy versions of a given set of reactants are omitted. Another group of 
possible reactions not included in table I are those resulting in three prod- 
ucts. The same reactants used above could also yield the following: 

CO + NO 5 CN + 20 917 kj/mole 

It is easily seen that the same arguments applied above may again be used in 
this case. Additional energy is required to produce the third product thereby 
allowing three-product exchange reactions to be ignored as possible controlling 
reactions . 

A third group of reactions which bear discussion are those containing spe- 
cies of primary interest but also requiring a species which is not significant 
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in determining the post-shock equilibrium state of the gas . An example is 
given by 


2CN £ N 2 + C 2 29 kJ/mole 

Even though this reaction requires relatively little energy, it necessitates 
admitting to the presence of C 2 . Since C 2 always appears in extremely small 
quantities under the circumstances of interest here, the CN lost or generated 
by this mechanism is trivial compared to other higher energy reactions using 
more predominate species. Thus, reactions involving such species as N^O, 
and NO 2 , to name a few, may be neglected and their kinetics studied using the 
thermodynamic profiles computed in their absence . In view of the preceding 
arguments, table I is believed to contain all the reactions of primary impor- 
tance to the chemical system for the shock-wave conditions of this study. It 
also contains a few which may be superfluous because of the nature of the 
system. These are identified in the "Results and Discussion" section by com- 
paring the relative effectiveness of all the reactions in producing a given 
species . 


Reaction Rates 

Before proceeding to an application of the postulated chemical system, it 
is first necessary to define the rate coefficients associated with each reac- 
tion. A review of recent literature (e.g., ref. 2) has shown that for air 
reactions of interest, a given rate is seldom known within a factor of 5 at 
high temperatures. Fortunately, there are some exceptions such as N 2 and 0 2 
dissociation rates which have been claimed accurate within a factor of 2 when 
collision partners are specified. However, in the cases where limited or no 
measured data are available, it can be considered exceptional to predict the 
correct rates within an order of magnitude. This "plus or minus one order-of- 
magnitude" criteria will be used to evaluate the theoretical predictions made 
in this study. The reactions given in table I may be separated for this 
purpose into two types, dissociation and exchange. Because of this separa- 
bility, it is advantageous to use the rate theory found most successful for 
each type . 

Dissociation reactions .- Because the dissociation reactions are assumed 
in this study to occur only with a general collision partner. It is necessary 
to use a rate theory that will aid in unifying the many published rates which 
specify collision partners. Also, a method is required to predict the rate of 
CN dissociation, the only dissociation reaction where no data were found. A 
recent and successful theory formulated by Hansen (ref. 12) from a version of 
the available energy theory is used for these purposes. Hansen 1 s results may 
be applied to the general dissociation reaction 


AB+M^ A+B+M (l) 

k f 
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The rate of production of molecule A, for example, by the above mechanism is 
defined as 

MAI = [M] J [AB ] _ A. [A][B]l (2) 

at o' l K c J 

where it is assumed that K c = kf/k^ 


o 


' 

1, A = B 

2, A ^ B 


and k^ = kf(T) 


The forward rate coefficient, kf, can he equated to that given hy Hansen. 
Rather than use Hansen* s results directly, however, two slight modifications 
were made to his expression for kf which allows it to he compared with the 
convenient Arrhenius form 


kf = PT^ exp(-E a /kT) 


(3) 


This is accomplished by neglecting the change in momentum cross section with 
increasing particle momentum and assuming D/kT » l/2. A difference of less 
than 5 percent is introduced in the final results at the temperatures of 
interest here. The modified Hansen equation for collision induced dissociation 
of a diatomic molecule, AB where AB ^ M, is then given hy 


k f 


'8k%^ 1/a ^Sqq A)\ 3/a 

- ^ ' 3^ce ' / 



w 


where e is the Napierian base. In equation (4), is equivalent to the 

Sutherland cross section used for viscosity (values given for several species 
in ref. 12) and D is the bond dissociation energy of molecule AB. Rewriting 
equation (4) using the mean values, = 30*10 _2Ora2 and p = l4, appropriate for 
diatomic molecules containing C, N, and/or 0, we have 


k f 


2X10 6 



D\ m 3 

kTy mole -sec 


(5) 


(It should he noted that the values of kf given in this paper are 10” 6 times 
those given in the more customary units of cm 3 /mole-sec . ) Equation (5) is 
compared with several published rate equations for 0 2 , NO, and CO dissocia- 
tion in figure 1. The published rates are plotted without giving consideration 
to their applicable temperature ranges although the composite experimental 
temperature limits of the individual experiments are indicated. Equation (5) is 
shown to provide rates which represent an approximate mean of the other values 
plotted; it therefore seems appropriate to use equation (5) for predicting the 
unknown CN dissociation rates. The difficulty in choosing among rate coeffi- 
cients where collision partners are specified is also avoided by using equa- 
tion (5) to define kf for all the dissociation reactions in table I except 
reaction 1 (C0 2 dissociation) . 
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Hansen (ref. 12) also derives a dissociation rate theory for polyatomic 
molecules such as C0 2 . When the results are compared with experimental data, 
however, they are not found to he so successful as for the diatomic case. 
Furthermore , several recent measurements of C0 2 dissociation rates given in 
references ik to 17 are shown by figure 2 to be in good agreement. Figure 2 
also includes Hansen 1 s unadjusted prediction which is shown to be high. The 
experimental results of Davies (refs. 16 and 17 ) were selected and have been 
used here as the C0 2 dissociation rate for reaction 1. 


Exchange react ions . - Of the eleven exchange reactions listed in table I 
only reactions 7 to 9 have rate coefficients based on experiment. Their 
values, taken from references 2 and 3; are given in table I. The remaining 
rates for the carbon reactions must be estimated. A review of the literature 
resulted in finding several methods for making such estimates but none showed 
substantial improvement over the classical quantum mechanical approximation 
described in reference l8. The method, commonly termed "Absolute Reaction Rate 
Theory" or "Trans it ion- State Theory," is also discussed with greater detail in 
the more recent references 19 and 20. The transit ion- state theory has an 
advantage over some other methods such as "collision theory" (see ref. 19 ) in 
that it contains greater molecular detail. In view of the structural similar- 
ity of the carbon, nitrogen and oxygen atoms which make up this kinetic system, 
the transition-state theory also lends itself to a convenient generalization 
based on a system for classifying the reactions, as will now be shown. The 
theory may be applied in one form to bimolecular reactions such as that given 
in general terms by 


AB + CD 



AC + BC 


( 6 ) 


where A, B, C, and D may be atoms or molecules. The asterisk refers to the 
unstable "activated complex" which serves as a transition molecule during the 
reaction. The forward reaction-rate coefficient, k^, is similar in definition 
to that used in equation ( 2 ) and is expressed by the transition-state theory 
as 


k« - 


kT £*_ 
h I1F V 


exp - 


kT 


(7) 


In equation ( 7 ), F is the partition function per unit volume of each reactant, 
r, and F* is thesame as F r but for the activated complex less the contribu- 
tion of one vibrational mode. The omitted vibration is the mechanism by which 
the complex is assumed to separate into reaction products. The transition- 
state theory does not provide a means of estimating the activation energy, E a . 
Such estimates will be discussed subsequently; here we are dealing only with 
the pre -exponential or "frequency factor" defined as 


A — 

h nF r 
r 


( 8 ) 
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The partition functions, F r and F^, are evaluated assuming equipartition of 
energies and contributions of molecular translation, rotation, and vibration 
only . Thus 


F 


fj^fy 

V 


( 9 ) 


where 


ft 


(^) 3/2 v 


(10a) 


f 


r 



(10b) 


V 



(10c) 


At this pointy one could substitute values of the individual molecular 
properties and evaluate A for each reaction in detail. However, this is not 
practical since these properties, which are well known for most reactants, are 
difficult to determine for an activated complex. Furthermore, equation ( 7 ) is 
derived by assuming a physical model of the reaction process which is highly 
idealized and also requires E a to be estimated. It therefore appears equally 
reliable to express A in a general sense for a given set of similar reactions 
(i.e., with similar reactants). The approach here will be to demonstrate the 
similarity of the C, 0, and N atoms and through this similarity construct a 
set of reaction-rate approximations based on the molecular configuration of 
pairs of reactants. To this end, table II lists values of m, 0 r , and Q v 
required by equations (10) for the species of interest. It can be seen that 
when properly grouped like quantities show a degree of similarity. For 
example, the monatomic species are of nearly the same atomic weight and, from 
the periodic table, have similar atomic structures. In the case of the 
diatomic species an average value of Q r or 6 V can be found which adequately 
approximates the thermodynamic properties of any of them. The same is true for 
the polyatomic molecules when the molecular arrangement is considered (i.e., 
linear or nonlinear) . The individuality of each species can therefore be 
dropped and attention given only to the number of atoms in each molecule and 
the manner in which they are arranged. This technique has been proposed in 
several references (e.g., ref. 21 ) but with different results since the con- 
stants chosen were for more arbitrary reactants. The results given here are 
specialized for C, 0, and N systems. 

Generalized rate coeff icients.- Consider the general bimolecular reaction 
given by equation ( 6 ). If equation ( 9 ) is substituted into equation ( 8 ) and 
the contribution from each reactant is denoted by an appropriate subscript 

a , = s — Y( ( tf y* (id 

eq.* ( 6 ) h (ft f r f v)A£( f t f r f v)cD 
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If now AB is an atom (called B) and CD a diatomic molecule , the reaction is 


B + CD -» 



products 


Since f r = f Vx2 = 1, equation (ll) becomes 


-v B 


- _ kT v ( f t f rfy)^. 

eq^.(aa) h f (f f f ) 

I t B ^ I t I r I v'CD 


( 12 ) 


( 13 ) 


Thus, A takes a special form for reactants of the molecular configuration in 
equation (12) . This scheme of reaction classification is now extended to other 
reactant configurations using suitable average values of the molecular 
constants . 


The translat ional energy contributions are treated first since all 
particles contain a translational energy mode. The translational terms , there- 
fore, may be factored from equation (ll) for any reaction. Incorporating 
equation (10a) equation (ll) is then rewritten as 


where the constant 


A = Cjl 


- 3 / 2 . 


T 


- 1/2 


(f r fv). 


^ f r f v^AB^ f r f v^CD 


C 


k 

h 



(l^a) 


(C = 6.68x 10 T for A in m 3 /mole-sec) and \i is a reduced molecular weight, 
an average molecular weight of l4 is used for all atomic particles, jl is 
approximated by 


11 = 14 


n AB n CD 
n A B + n CD 


If 


(l4b) 


The rotational partition functions are approximated using equation (lOb) with, 
average values of Q . From table II the average values selected are 


Diatomic molecules 6 , 


3 - 3 ° K 


Polyatomic molecules Q m 0.6° K 

When the proper values of l are substituted into equation (lOb) 

1 Atoms l = 0 

( 0.3T Diatomic and linear molecules l = 2 


2T 372 Nonlinear molecules 


l = 3 


( 15 ) 
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If 0v a / T is assumed small , the vibrational contributions as represented by 
equation ( 10 c) can be approximated by 


L v 


s 

n 


a=i 


( 16 ) 


a 


The values of Q Va given by table II are replaced by average values selected 


as 


e Vl = 2900 ° k 

CD 

0v 2 = 1900° K 

0 v 5 

e Y = 930° k 

v 3 

0 v a 


0 V = -4*300° K (930° K for linear molecules) 
0 V5 = 4300° K 


The number of vibrational degrees of freedom, s, must still be determined. 

The value assigned to a particular molecule depends on the number of atoms , n, 
in the molecule. Since an activated complex has one less vibrational mode 
than its molecular counterpart, the following rules apply: 

3n-5 Linear or diatomic molecules 

s = { 3n-6 Linear complex or nonlinear molecule 

3n-7 Nonlinear complex 

All activated complexes are assumed nonlinear and the vibrational mode corre- 
sponding to 0 Vl is considered the mode causing destruction of the complex. 
From the preceding, f v is then given for the configurations needed as 


f V =< 


r 3*4xiO“ 4 T 

Diatomic molecule 

5.7xict 7 't 2 

Triatomic complex 

2. 0X10“ ^T 3 

Nonlinear triatomic molecule 

2 .1X10“ ^T 4 

Linear triatomic molecule 

6 . 1 x 10 “ ^T 5 

Four atom complex 

4 . 9xio - 29 T 8 

Five atom complex 


( 17 ) 


The combined results of equations (14) to ( 17 ) are then as follows: 
Denoting an atom as A, a diatomic molecule as A^B or C'D, a linear triatomic 
molecule as C*D*E, and a nonlinear triatomic molecule as D^E, we have 
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Reactant pair 


A ( m 3 /mole - se c ) 

A+B-C 2X10 4 T 

A+C -D -E 3X10 2 T 

A+D C E 5X10 _2 T 3/2 

A-B+C-D 1X10-3T 2 

A-B+C-D-E 2X10- 5 0 P 2 

A-B+D C E 3X10 _9 T 5/2 

These are the values of A used for reactions 10 to 17 in table I. 

Comparisons with experiment .- Figure 3 provides comparisons of equa- 
tion (l8) with values of A obtained by experiment or other estimates. The 
applicability of figure 3 relies on the assumption that the Arrhenius form of 
kf (eq. (3)) is valid over any temperature range. It is generally accepted 
that this is not the case but since all of the experimental results plotted in 
figure 3 were obtained by fitting equation ( 3 ) to the basic data, figure 3 
should at least provide an indication of the theoretical accuracy. The^figures 
are therefore intended as a comparison of the predicted magnitudes of A with 
experiment in the experimental temperature range (log T ^ 3) and also as a 
comparison of the functional dependence of A on temperature as given by 
theory with that indicated by experiment (i.e., by comparing the slopes of the 
lines plotted) . The experimental values plotted are extrapolated considerably 
beyond the temperature limits of their associated experiment. The maximum 
temperature of most experiments represented did not exceed 5000° K. The three 
general reactions shown are those used for this study. The shaded box on each 
figure indicates the range of nonequilibrium temperatures covered for 
V s = 8 km/sec and the assumed inaccuracy of A. The specific comparative rates 
are from the references given in parentheses beside the reaction. Figure 3 
relies on the premise that a method successful in predicting N, 0 reactions 
works equally well for C, N, 0 reactions because of the similarity of the 
carbon atom to the nitrogen and oxygen atoms. This premise is necessary since 
only N, 0 reaction rates were available for comparison. Figures 3(a) and 3(b) 
show excellent agreement between theory and experiment considering the gross 
assumptions made. A probable error factor of plus or minus an order of mag- 
nitude is assigned for further analysis of the reaction- rate uncertainties. 

This error factor is based on the experimental results shown although larger 
errors may exist. In figure 3(c) the scatter in experimental values does not 
allow any definite conclusions regarding the accuracy of the theory. The 
experimental values for a single reaction are not in good agreement (viz., 

^ 2+02 2N0) . An error factor of plus or minus two orders of magnitude is 

therefore assigned to reactions of this type. Fortunately, such reactions 
involving two diatomic reactants can be shown to be minor contributors to this 
kinetic system unless their rates exceed the present estimates by four or five 
orders of magnitude. It is seen from figure 3(c) that this is conceivable but 
it is believed “unlikely. 
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Activation energy An estimation of k^ using equation ( 7 ) still 
requires a value of the activation energy., E a . A method of estimating E a 
using intermole cular potential energy surfaces is described in reference l8 
but it requires elaborate computing techniques which are not warranted con- 
sidering the approximate nature of this study • Some semiempirical rules have 
also been proposed by Glasstone, Laidler, and Eyring (ref. l8) and Hirchfelder 
(ref. 22) but were adjusted for reactions involving hydrogen and do not give 
good results for the N, 0 reactions. An empirical scheme was therefore con- 
structed which follows the basic ideas of reference l8 but is centered on N, 0 
reactions (assuming again that C, N, 0, and N, 0 reactions are similar). This 
scheme results in an empirical relation which is general for all types of 
bimolecular reactions used here and is written for an endothermic process. 

To facilitate a discussion of types of reactions, the following endo- 
thermic reaction types are identified: 


1 

AB-HVl -*• A+B+M 

dissociat ion 

II 

AB+C A+BC 

three center exchange 

III 

AB+CD -* AC+BD 

four center exchange 


where each symbol is an atom or pair of atoms not separated by the reaction. 

The following postulates are proposed: 

(1) For Type I reactions, E a = Hr = D^-g. validity of this has been 

demonstrated experimentally with a few exceptions such as C0 2 dissociation. 

It is also reflected in equation (k) for dissociation reaction-rate 
coefficients . 

(2) For any reaction occurring directly from ground-state reactants, 

E a > Hr. Measurement verifies this for all of the exchange reactions examined, 
and it is a requirement from an energy conservation viewpoint. 

(3) E a is a function of the sum of the bond dissociation energies for 
those bonds which are broken in the reaction. This presumption is a result of 
the methods suggested in references l8 and 22. It may be given a pseudo- 
physical basis by relating it to the energy required to destroy the activated 
complex and neglecting its heat of formation. 

Regardless of the lack of rigor in the preceding postulates, they appear 
to represent the behavior of experimental observations and are therefore 
pursued. Normalizing with respect to Hr, the third postulate gives 



The boundary conditions become E a /Hp = 1 when (ZRr)/H-R = 1 and E a /HR > 1 
always. Equation ( 19 ) suggests a correlation of experimental results as given 
by figure k. The data are taken from reference 2 and pertain to N, 0 reactions 
except where noted otherwise on the figure . A linear correlation is obtained 
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on semilog axes if the likely degree of error in the experimental data is kept 
in mind and temperature dependencies are ignored. From figure b, the following 
empirical relation is then proposed: 


- exp 

% 


1 f £Dr 

10 W 



( 20 ) 


The error assigned to the results of figure b requires a multiplying factor for 
the right-hand side of equation (20) of exp ± 3/l0. For values of 
H R /k < 50,000° K, this error affects k f hy less than an order of magnitude 
(Hft/k is less than 50 , 000 ° K for most of the carhon exchange reactions and all 
those of importance) . Equation (20 ) was used to complete the estimates of kf 
for reactions 10 to 17 of table I. 


Computing Method 

The nonequilibrium properties behind a normal shock wave were computed 
using a modified version of the IBM 70^ program described by references 23 and 
2b . The system of equations will not be discussed in detail here because of 
the completeness in description supplied by reference 2b . However, in brief, 
the program relies on conventional gasdynamic and chemical conservation equa- 
tions in conjunction with a quantum-mechanical description of the individual 
species thermodynamics. Thus, spectroscopic data such as that given in 
table II are used as input . Some additional input data also required were the 
heats of formation of the individual species and were obtained from refer- 
ence 25* The computer program given by reference 2b was recoded for use on an 
IBM 709^ and modified to accept polyatomic molecules with up to nine vibra- 
tional degrees of freedom. Equilibrium electron excitation has also been 
included as in the original program. Electronic energy levels up to 
80,000 cm™ 1 (l. 6X10" 18 J) were incorporated although such refinement has only 
small effects on the numerical results of this study. 

Computing time was found to have a much stronger dependence on the number 
of species than on the number of reactions. Minor species retained in the 
calculations slow the computing progress appreciably and, in some cases, their 
computed concentrations suffer from excessive numerical errors. Typical com- 
puting times for a relaxation to within 10 percent of equilibrium at 
V s = 8 km/sec were about 5 minutes. 

RESULTS AMD DISCUSSION 


With all the necessary reactions and rate coefficients established, this 
section is devoted to an analysis of the chemical kinetics for shock condi- 
tions described previously. A discussion is first given of some results using 
the nominal reaction rates in table I. General characteristics of the chemical 
kinetics are illustrated and several features which are independent of the 
reaction-rate errors are discussed. This is followed by an evaluation of the 
changes in thermo chemical profiles due to uncertainties in the reaction rates. 
The consequences of deleting certain species and reactions from the chemical 
system are also investigated in an effort to simplify the kinetic model. 
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Finally a brief study is made of the CN kinetics over a range of shock speeds 
because of the particular interest in that species for studies of shock-layer 
radiat ion. 


General Characteristics 

Some estimates of the nonequilibrium chemical behavior for a shock speed 
of 8 km/sec using nominal reaction rates are presented in figures 5 to 7. 
Thermodynamic profiles and species concentrations are plotted as functions of 
the distance parameter y, which is equivalent to a binary scaling parameter 
as described by reference 26. The parameter y may be converted to laboratory 
time by the approximation 

t „ (21) 

L (pi/po)V s 

which is correct within 10 percent for all shock conditions considered. 

Part (a) of each figure illustrates profiles of the thermodynamic quantities 
most sensitive to the chemical kinetics. Other thermodynamic variables which 
are not shown, such as local pressure or enthalpy, vary only a few percent 
behind the shock wave. Local velocity is also omitted since, for a normal shock 
wave, it is easily obtained from 


v = 4 (22) 

P 

Part (b) of each figure illustrates the individual species concentrations 
given as mole fractions of the ambient gas mixture. These may be converted to 
the more familiar mole fractions based on local gas composition by 

E. = my. (23) 


The equilibrium thermodynamic and species concentration values given in fig- 
ures 5 to 7 were obtained using a computer program developed by Simcox and 
Peterson (ref. 9 )* The equilibrium species concentration levels indicated 
separately on the figures are for P±/p Q ~ 10 -6 only since the nonequilibrium 
curves for p-J p Q = 10” 1 achieve equilibrium levels within the scale of the 
figures. In some cases, slight differences between the fundamental inputs 
used in reference 9 and in this study preclude a close comparison of equilib- 
rium species concentrations with the nonequilibrium profiles. However, the 
discrepancy never exceeds 5 percent. Tabulated equilibrium concentration 
values are given in table III to provide those numerically too low to appear 
on the figures. 

Figures 5 to 7 demonstrate several features of the chemical kinetics 
which are retained regardless of the reaction rates used. For example, the 
profiles, when plotted in terms of y, show no dependence on ambient density 
until equilibrium is nearly achieved. This occurs because the dissociation 
reactions proceed predominantly in the forward direction (except near equilib- 
rium) while the exchange reactions are binary in either direction. Among other 
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things, this correlation allows easy estimation of the time or distance to 
equilibrium for intermediate ambient densities. That is, if the equilibrium 
values of any quantity sensitive to the kinetics are known, the manner in which 
its values depart from the lowest density curve to approach equilibrium is not 
difficult to estimate. (Equilibrium thermodynamic properties and species con- 
centrations for initial gas mixtures close to those used here are given by 
ref. 11.) A second feature demonstrated in figures 5 to 7 pertains to the 
manner in which those species quickly generated (e.g., NO and CN) tend to 
overshoot their equilibrium concentrations . A “peak” concentration for each 
generated species may be defined whether it is actually achieved or not. It 
depends only on the shock speed and ambient gas composition, assuming reaction 
rates are known. This density- independent peak is achieved if the equilibrium 
concentration is less than the peak value. In the event that the equilibrium 
concentration is greater, the concentration will rise monotonically . An 
example of both cases can be found in figure 6(b) by comparing the CN concen- 
trations for p / p Q = 10 -6 and 1CT 1 . Notice in the case where p x /p 0 = 10 -6 , 
the CN concentration reaches a peak and then decreases to its equilibrium 
value. This same behavior would occur for all lower densities and any higher 
density where the equilibrium CN concentration was less than the "peak." For 
those ambient densities where the equilibrium CN concentration was higher than 
the "peak," the CN concentration profile would rise monotonically as it does 
for Pi/ P 0 = 10” i shown in figure 6(b). By virtue of this same behavior, 
certain minor species which are strong radiant emitters (e.g., C 2 ) were also 
found to achieve concentrations considerably higher than their equilibrium 
value within the nonequilibrium region. 

A comparison of the (a) parts of figures 5 to 7 to obtain the effects of 
ambient gas mixtures showed only minor differences for the mixtures considered. 
The differences in thermodynamic profiles arise principally because of small 
changes in their end points (i.e., the initial or frozen and equilibrium 
values) . The time or distance required to reach near-equilibrium is essen- 
tially unaltered. Conversely, comparing part (b) of each figure shows ambient 
gas composition to have sizable effects on the species concentration but, again, 
the major differences are in the relative initial and equilibrium values rather 
than time or distance. However, the example is cited where decreasing the 
ambient gas composition from 50 to 5 percent C0 2 causes approximately a four- 
fold increase in y at the CN peak. On the other hand, y at the NO peak is 
unaffected. 

To gain insight into the relative importance of each reaction to the whole 
chemical system it is worthwhile to compare the rate of production of certain 
species by individual reactions. Figure 8 makes these comparisons at one set 
of shock- wave conditions for all except the monatomic species. The dependent 
parameter, Q. ., in figure 8 is proportional to the species production rate 
defined by equation (2); that is, Qij in figure 8(a) represents the net pro- 
duction (or reduction) of CO by each reaction shown and accounts for both the 
forward and reverse contribution of that reaction. Some reactions (viz., 8 and 
13, table i) show a negligible contribution to any total species production. 
Reaction 13 would not become significant unless kf for that reaction was 
several orders of magnitude greater than those values regarded as reasonable 
(i.e., using the limits defined by fig. 3(c))* Reaction 8 has a relatively 
well-established rate which keeps it ineffective for any species production, 
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thus making it clearly unimportant (e.g., figs. 8(b) and 8(f)). Some other 
reactions appearing in figure 8 also show the possibility of being deleted. 

For example, reaction 16 (figs. 8(a) and 8(e)) is minor compared to other 
reactions contributing to the same species. Also, it is a reaction most 
probably already contained in the measured C0 2 dissociation rate (reaction l) 
since it incorporates only products of C0 2 . Its effect on the kinetics has 
been numerically tested and found to be negligible. Reactions 15 and 17 
(figs. 8(a) and 8(f)) also incorporate only C0 2 and its products. However, 
they serve as principal 0 2 generators and cannot be deleted when considering 
that species. 

Figure 8 is representative of similar plots for other ambient densities 
and gas mixtures. The results of figure 8 are altered by higher ambient 
densities in a manner compatible with that of figures 5 to 7; that is, Qj_j 
for a higher density departs from the low density curves towards zero when 
equilibrium is approached. _ A change in gas mixture will cause some slight 
relative shifting of the Qj_j curves but does not alter the conclusions 
obtained regarding the importance of individual reactions. Hence, reactions 8, 
13 ; and 16 may be deleted for all ambient densities, gas mixtures, and, in 
addition, other shock speeds. 

Additional thermodynamic and species concentration profiles for V s = 6 
and 4 km/sec are presented in figures 9 and 10, respectively. These figures 
illustrate the diminishing differences between the frozen and equilibrium 
values of any quantity with decreasing shock velocity and show a lessening 
tendency of the generated species to reach high peak concentrations. The 
extent to which the nonequilibrium region is lengthened with decreasing shock 
velocity is also demonstrated. The effects of other ambient densities and gas 
compositions at these shock speeds are similar to those at V s = 8 km/sec. 


The Effects of Reaction-Rate Uncertainties 

When the estimated errors which may be present in rate coefficients for 
the carbon exchange reactions are recalled (figs. 3 and 4), it is of interest 
to investigate their effects on the chemical behavior as given in figures 5 to 
7* For this purpose, the shock conditions of figure 7 were chosen as an 
example since they are the most sensitive to reaction-rate uncertainties of all 
the cases shown. From the limits given by figures 3 and 4, rate combinations 
were chosen which would maximize or minimize the generated diatomic species 
concentrations since their behavior is the most sensitive to exchange reactions. 
To obtain a "maximum profile,” maximum values were used for the rate coeffi- 
cients of all reactions in question except 10. Figure 8(c) shows reaction 10 
quickly reverses, thus it must be a minimum for maximum CN formation. The 
"minimum profiles" came from opposite considerations. 

The resulting thermodynamic profiles, given by figure ll(a), are rela- 
tively insensitive to the rates used. In particular, local density is the 
variable most insensitive making its measurement inadequate for resolving 
chemical detail in a complicated system such as this one (density is commonly 
measured by interferometer or electron beam techniques to obtain rate informa- 
tion) . In contrast to the thermodynamic quantities, figures 11 (b) and ll(c) 


18 



show species concentrations to exhibit high sensitivity to the rate variations. 
This is particularly evident where the time or distance to a peak concentration 
and its associated magnitude are examined. Thus y in the absence of an ability 
to resolve chemical detail by other methods., considerable information may be 
gained by spectroscopically measuring the time to peak CN concentration for 
example. This technique would not provide individual reaction rates nor define 
the underlying mechanisms in the customary detail but it would allow a more 
refined reaction system model to be established; suitable for the needs of many 
predictions in aerodynamic applications . 


Contributions of the 0 2 Reactions 

The analysis of experimental data for a chemical system of many species 
and reactions is always made simpler by removing as many insignificant varia- 
bles as possible from the theoretical model. When only the gross thermo- 
chemical behavior is of interest; removal of any species or reaction reduces 
the number of reaction rates which must be measured and aids in reducing the 
numerical difficulties encountered when calculating such behavior. 

The significance of the contributions of 0 2 reactions to the remaining 
kinetic system becomes questionable when the results in figure 8 are examined. 
As seen there; none of the 0 2 reactions are important to the formation of any 
species except 0 2 . Thus; the five reactions containing 0 2 were deleted and the 
results compared to nonequilibrium profiles including them. Figure 12 gives 
one example of such a comparison. The case shown (50 percent C0 2 ) is that dis- 
playing the greatest effect of all the ambient mixtures considered using nomi- 
nal reaction rates. The thermodynamic properties (fig. 12(a)) are shown to 
have a negligible dependence on 0 2 while the effect on concentrations of other 
species is small (fig. 12(b)). These small effects are not enlarged when other 
reaction rates are used except in the extreme case where the 0 2 reaction rates 
are maximized and the others put to a minimum. However; that situation is not 
seriously considered since experimental evidence favors the opposite case. 2 It 
is interesting to note that Howe and Sheaffer (ref. 4) also found the inclusion 
of 0 2 unnecessary when considering 100 percent C0 2 . 

The effect of deleting the 0 2 reactions in a 50 percent C0 2 mixture was 
also investigated at V s = 4 km/sec where 0 2 appears in larger equilibrium 
concentrations. The effects at that speed are larger but; again; they appear 
primarily in the species concentrations. Thus ; for shock speeds greater than 
4 km/sec and reaction rates equal to or greater than those used here; 0 2 may be 
omitted without serious error. The results of such a deletion are a 20- to 
40-percent reduction in computing time and a simplification of the reaction 
system from l4 to 9 reactions; assuming reactions 8; 13; and 16 have already 
been removed. 


^Experimental evidence” refers to some unpublished preliminary measure- 
ments of CN radiation intensities made by James Arnold of the Ames Research 
Center. The data were taken behind incident shock waves in a shock tube at 
conditions close to those considered here and indicated a relatively rapid 
formation of CN . 
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Cyanogen Kinetics 


Special attention is given here to the kinetics of CN because of its 
importance to the total radiant emission from behind shock waves in these gas 
mixtures. Figure 13 gives an estimate of the CN concentration at its peak 
value and the laboratory time to reach the peak for a range of shock speeds and 
two gas mixtures. The nominal reaction rates were used for this calculation 
and the 5 percent C0 2 mixture was omitted because of its close proximity to the 
10 percent C0 2 case . The possible error in this estimate may be evaluated by 
considering the variations in definition of the CN peak as shown in fig- 
ure ll(b) . The percentage of shifting due to reaction-rate uncertainties shown 
in figure 11 (b) at V s = 8 km/sec is essentially preserved at the lower speeds. 
The effect of ambient gas composition on CN peak concentration can be seen from 
figure 13 to be fairly significant. Figure l4 compares the peak (as given in 
fig. 13) with the equilibrium CN concentrations for the densities used in this 
study. It is of interest to note (cf., figs. 14(a) and 14(b)) that for a 
given shock speed the 50 percent C0 2 mixture has a greater tendency to achieve 
a peak CN concentration than the 10 percent C0 2 mixture; that is, the magni- 
tudes of the peak concentrations lie above the equilibrium concentrations for 
a greater range of ambient densities in the 50 percent C0 2 case. Thus, non- 
equilibrium CN radiation, for example, at a given ambient density and shock 
speed, may be more significant for the higher initial C0 2 concentration mix- 
tures at least up to 50 percent C0 2 and particularly at shock speeds less than 
8 km/sec. This may be compared to conclusions drawn from equilibrium results 
(see, e.g., ref. 7) which show maximum CN equilibrium concentrations to occur 
when the initial C0 2 content is 9 to 16 percent. 

Finally, attention is called to the trends shown in figure l4(b) at 
V s < 5 km/sec. The extreme overshoot of the CN concentration at these lower 
speeds is typical in regions where the amount of CN present is most sensitive 
to temperature. It can be seen from the equilibrium curves that small 
increases in V s (or correspondingly, in temperature), in the region of 
V s = 5 km/sec, produces large increases in equilibrium concentration. Thus 
there is, in effect, a "cutoff temperature" below which the CN concentrations 
are negligible. If, as in the case for V g < 5 km/sec, the equilibrium temper- 
ature is less than the cutoff temperature but the nonequilibrium temperatures 
are greater, there will be extreme overshooting of the CN concentration in the 
nonequilibrium region. However, it should also be noted that although the 
ratio of peak to equilibrium concentrations is large at these conditions, the 
magnitude of the peak concentration is still much lower than those achieved at 
higher shock speeds. Thus, the effects of nonequilibrium phenomena on CN radi- 
ation intensities (for example) may be more pronounced at shock speeds less 
than 5 km/sec and at a given density but the magnitude of the intensity will 
also be very low. Similar behavior at higher shock speeds is found for other 
minor species such as C 2 , N2O, and N0 2 . 


CONCLUDING REMARKS 


A study of the chemical kinetics behind normal shock waves in mixtures of 
C0 2 and N 2 has been made as a prelude to experimental investigations. Emphasis 
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was placed on shock conditions typical of those anticipated during entry into 
the Martian atmosphere. A shock speed of 8 km/sec was studied in detail hut 
lower speeds are also discussed. Ambient gas mixtures containing 10, and 
50 percent C0 2 and ambient densities of 10 -6 and 10~i times that of sea level 
were considered. 

A survey of current methods of estimating reaction rates became necessary. 
It was found that most of the dissociation reaction rates could be based on 
measured data and that an adequate theory was available for estimating those 
that could not. Most of the exchange reactions required estimated rates and 
the "transition-state" approximation was used. The results were compared with 
measured data where available, and agreement was found to be satisfactory. A 
scheme for correlating activation energies and predicting unknown values was 
also devised for reactions containing C, N, and 0 atoms. 

The species considered were C0 2 , N 2 , CO, NO, CN, 0 2 , C, 0, and N. An 
original system of 17 bimolecular reactions was reduced to ik- for the species 
above. Deletion of 0 2 and its reactions introduced only small errors and 
further reduced the chemical system to 9 reactions. They are 

C0 2 + M 5 CO + 0 + M 
N 2 + M X 2N + M 
WO+M^N+O+M 
C0+M^C+0+M 
CN+M^C+N+M 
N 2 + 0 5 NO + N 
N 2 + C ^ CN + N 
CO + N X CN + 0 
CO + NO X C0 2 + N 

The remaining 0 2 reactions, important mainly to the production of 0 2 
itself are 


0 2 + M X 20 + M 
NO + 0 ^ 0 2 + N 
C0 2 + 0 CO + 0 2 
CO + 0 t o 2 + c 
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The chemical kinetics were found to allow the assumptions of binary 
scaling over most of the nonequilibrium region. This led to a correlation of 
the dependence on ambient density and allowed density- independent results to be 
presented except near equilibrium. 

Certain generated species such as CO, NO, CN, and 0 2 were shown to over- 
shoot their equilibrium concentrations for most ambient densities of interest 
in the flight regime. This was particularly the case for CN which suggests 
significant nonequilibrium radiant emission at low densities. 

The effects of ambient gas composition on the length or duration of the 
nonequilibrium region behind the shock wave were shown to be small. The prin- 
cipal differences in thermochemical profiles between the several ambient gas 
compositions are due to changes in relative magnitudes of the profile end 
points rather than the time or distance required to reach near-equilibrium. 

Admitting errors of plus or minus an order of magnitude in the estimated 
reaction rates caused only small differences in the predicted thermodynamic 
profiles. The species concentration profiles, however, are very sensitive to 
such error. Measurement of the time to peak concentration of an identifiable 
species was shown to be a possible means of establishing rate information. 

Special attention was given to the kinetics of CN. It appears to be one 
of the species most sensitive to ambient gas composition and reaction-rate 
uncertainties . CN displayed a greater tendency to overshoot its equilibrium 
concentration in increasing ambient proportions of C0 2 to 50 percent C0 2 • 

At shock speeds less than 5 km/sec, the CN concentrations are small in the 
nonequilibrium region but still several orders of magnitude greater than their 
equilibrium value. The extent of concentration overshoot decreases with 
increasing shock velocity but remains significant at the shock speeds of 
interest for entry into the Martian atmosphere. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field , Calif., Oct. 8, 1965 
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TABLE I.- PRELIMINARY REACTION SYSTEM 






Assumed reaction- 





rate terms* 

No. 

Group 

Reaction 

A 

% (°K) 

1 


k 

C0 2 +M 5 C0+0+M 

1.2X10 5 T 1/2 

3^360 

2 

1 

c 

4 

4 

) 

N 2 +M 5 2N+M 

8 . 6X10 ^T -1 

113350 

3 

•r 

+■ 

a 

3 

-t 

NO+M 5 N-K)4M 

2 ^XIO 13 ! -1 

75560 

4 

CJ 

O 

CQ 

rn 

C04M ^ C+O+M 

S^XIO 13 ! -1 

I 2896 O 

5 

•r 

f= 


CN+M ^ C+N4M 

5.2X10 13 T _;l 

94i4o 

6 

2 

r 

0 2 +M 5 20+M 

2 ^XIO 13 ! -1 

59390 


A 

K 




7 



N0+0 5 0 2 +N 

3 .2X1 0 3 T 

19680 


U 

§ 




8 

•H 

x 

o 

X 

Cl) 

N 2 +0 2 5 2N0 

2 .0X10 8 

61600 

9 


f 

N 2 +0 5 NO+N 

7-0X10 7 

38000 

10 

1 

i 

N 2 -tC CN+N 

2X10 4 ?1 T 

31560 

11 



CO+N ^ CN+0 

2X10 44l T 

45800 

12 



1 

co 2 +n ^ cn+o 2 

3X10 2:i:;i T 

49560 

13 

a 

o 

<L> 

N 2 +C0 5 CN+NO 

lXlO^+ST 2 

92010 


X 

eg 




l4 

u 

ctf 

o 

X 

CJ 

X 

<D 

CO+NO 5 C0 2 +N 

1X10 _ 342 T 2 

20980 

15 


C0 2 +0 ^ C0+0 2 

3X10 2:F1 T 

18210 

16 



2C0 5 C0 2 +C 

lXlO -3 * 2 ! 2 

72390 

17 


r 

C0+0 5 0 2 +C 

2X10 4 + 1 T 

69570 


*k = A exp (- ^ m3 
f kTy mole-sec 
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TABLE II.- SOME MOLECULAR CONSTANTS 



CD 

Species 

A 

m 

e r (°K) 

a=l 

a =2 

a =3 

a =4 

c 

12.011 






0 

16.000 






N 

l 4 .007 






CO 

28.011 

2.78 

3122 




NO 

30.007 

2.45 

2740 




CN 

26.018 

73 

2976 




n 2 

28.014 

5.78 

3395 




0 2 

32.000 

4.16 

2274 




c 2 

24.022 

4.70 

2362 




*co 2 

44.011 

1.12 

3380 

1928 

960 

960 

*n 2 o 

44.014 

0.60 

3199 

1849 

847 

847 

**no 2 

46.007 

0.36 

2332 

1899 

932 



*Linear molecule 
**Nonl inear molecule 
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TABLE III.- EQUILIBRIUM SPECIES CONCENTRATIONS (7.) FOR V g = 8 km/sec 

J 


i Ambient 
mixture 

5 percent C0 2 ■ 

- 95 percent N2 

• 10 percent C0 2 

90 percent N 2 

i 

50 percent C0 2 

- 50 percent N 2 

p i 

Po 

10 _s 

H 

1 

O 

H 

10~ 6 

H 

1 

O 

H 

CO 

1 

0 

H 

10-1 

C 


3.60-2 

7.87-2 

6.49-2 

2.35" 1 

I.6O-1 

0 

9.13- 2 

9.06-2 

1.77-1 

1.74-1 

7.35- 1 

6.69- 1 

N 

1-35 

9-46" 1 

1.26 

8.64-1 

7.99-1 

4.90-1 

CO 

5 .6l -3 

4.75" 3 

2 .17-2 

1.96-2 

2.64-1 

3-03” 1 

, NO 

9-6 7" 5 

4.34-2 

i.67“ 4 

8.11-3 

3.50- 4 

2 .26-2 

CN ' 

8.05- 4 

9.78-3 

1.48-3 

1.84-2 

2 .62-3 

3. 69- 2 

0 2 

| 

3.90- 8 

1 1 

I.94- 5 

1.50-7 

7.74“ 5 

2 .6l -6 

1.28-3 

n 2 

2.80 -1 

H 

o\ 

[>* 

H 

1 

OJ 

4.64-1 

1.02-1 

2.32-1 

co 2 

2 .12-9 

3.20-7 

1.03- 8 

2 .6l- 6 

2.67 -7 

I.86- 4 

m 

0.564 

0.639 

0.551 

0.623 

0.469 

0.522 


Note: A number such as 0.0025 appears as 2 -5 3 in the table. 
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Ref 3 , composite 



(c) NO dissociation. 


Ref. 13, experimental 
temperature range 



(d) CO dissociation. 


Figure 1.- Concluded. 




log, 0 A log lo A 


Transition state theory 
Measured values or 
other estimates 
Region used in this study 
for V s = 8 km/sec 



4 L 


(a) Reactions of the type A+B-C products. 

10 - 



(b) Reactions of the type A+B-C-D products. 

14- 


N 2 + 0 2 -2N0 



Io9io t 

(c) Reactions of the type A-B-tC-D -*■ products. 


Figure 3.- Frequency factors for some exchange reactions. 








<J0 

VJI 


(b) Species concentrations. 
Figure Concluded. 






(a) Thermodynamic properties. 

Figure 7 .- Nonequilibrium profiles; V s = 8 km/sec in 50 percent C0 2 - 50 percent N 2 









\ 



N 2 producing reaction 

N 2 using reaction 

Reactions too low to be shown: 
(8) N 2 + 0 2 = 2NO 
© N 2 + CO = CN + NO 


\ 

\ 

\ 



(d) N 2 production. 
Figure 8.- Continued. 
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I0 10 -' 


<D co 2 +m==co+o+m 


C0 2 producing reaction 

— C0 2 using reaction 


© C0 + N0=«C0 2 + N 
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I0 8 - 


I0 7 L 
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(e) C0 2 production. 
Figure 8 .- Continued. 









(a) Thermodynamic properties. 

Figure 11.- The effects of uncertain reaction rates; V s = 8 km/sec in 50 percent C0 2 - 50 percent N 








(a) Thermodynamic properties. 

Figure 12.- The effect of deleting 0 2 reactions; V s = 8 km/sec in 50 percent C0 2 - 50 percent N 2 

Pj/pQ = IQ"®. 




! 0‘ 7 _ 10 ' 
y 


(b) Species concentrations 


Figure 12 . - Concluded . 







V s , km/sec 


(b) 50 percent C0 2 - 50 percent N 2 
Figure l4-.- Concluded. 
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